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Abstract 

The time-dependent Ginzburg-Landau (TDGL) equation for a single component non- 
conservative structural order parameter is used to study the spatio-temporal evolution of a second 
phase in the vicinity of an edge dislocation in an elastic crystalline solid. A symmetric Landau 
potential of sixth-order is employed. Dislocation field and elasticity modify the second-order and 
fourth-order coefficients of the Landau polynomial, respectively, where the former makes the co- 
efficient singular at the origin. The TDGL equation is solved numerically using a finite volume 
method, where a wide range of parameter sets is explored. Computations are made for tempera- 
tures both above and below the transition temperature of a defect-free crystal T c q. In both cases, 
the effects of the elastic properties of the solid and the strength of interaction between the order 
parameter and the elastic displacement field are examined. If the system is quenched below T c o, 
steady state is first reached on the compressive side of the dislocation. On the tensile side, the 
growth is held back. The effect of thermal noise term in the TDGL equation is studied. We 
find that if the dislocation is introduced above T c q, thermal noise supports the nucleation of the 
second phase, and steady state will be attained earlier than if the thermal noise were absent. For 
a dislocation-free solid, we have compared our numerical computations for a mean-field (spatially 
averaged) order parameter versus time with the late time growth of the ensemble-averaged order 
parameter, calculated analytically, and find that both results follow late time logistic curves. 

Keywords: Heterogeneous nucleation; phase ordering; edge dislocation; Ginzburg-Landau the- 
ory; phase-field method 
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I. INTRODUCTION 



Defect enhanced nucleation of a new phase has been observed and modeled in a variety 

no 

of crystalline materials through the years 11H7J. In metallurgical systems, for example in 
aluminium-zinc-magnesium alloys, dislocations induce and enhance nucleation and growth 



of the coherent Laves phase MgZn 2 8Ml 01] . In ammonium bromide (NH 4 Br) a new phase 
appears in the proximity of crystal dislocations In titanium and zirconium alloys, 

containing hydrogen, when the hydrogen content exceeds the terminal solid solubility at a 
given temperature, hydrides (TiH x , ZrH^) form near dislocations and on grain boundaries, 



affecting mechanical properties of the alloys 12l.ll3j|. In magnetic systems such as gadolinium 



iron garnet GdsFesOi^, dislocations not only affect the spin reorientation transition lines in 
the associating phase diagram, they also change the nature of phase transition 

HQ. 

Other examples on the effect of dislocations include quantum crystals, metals and semi- 
conductors, in which electrons get localized near an edge dislocation producing discrete 
levels in certain part of the energy spectrum locally, hence altering the electronic structure 
of the solid, thereby affecting its thermodynamics and transport properties {ltsl \v\ . In a 



superconductor and a Bose solid, Ginzburg- Landau type model 



calculations show that 



superconductive and superfluid phase nucleate first (at higher temperatures) on dislocations 



prior to the bulk of the crystal 
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2l|. 



A generic theoretical approach to model the aforementioned phenomena is the Landau- 
Ginzburg theory of phase ordering. In this paper we use this theory to describe the kinetics 
of second-phase nucleation in the vicinity of an edge dislocation in solid crystals. A scalar 
non-conserved order parameter characterizing the presence or the absence of the new nucleus 
is assumed. The elastic behavior of the solid is taken into account by including the so-called 
striction term in the system free energy, which considers the interaction between the order 
parameter and deformation 
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23] . Moreover, the effect of a dislocation is incorporated in 



the system free energy through the equation of mechanical equilibrium in the presence of an 
isolated edge dislocation. The time-dependent Ginzburg-Landau (TDGL) equation with a 
symmetric tricritical potential energy is solved numerically to evaluate the spatial-temporal 
behavior of the order parameter in the vicinity of an edge dislocation for different sets of 
the coefficients in the potential energy. Our model is applicable to systems where second- 
phase ordering occurs with a preferred orientation of nuclei under an external force. These 



comprise a"-phase formation in Fe-N alloys 24J, #'-phase nuc 



hydride formation in titanium 



271 ] and zirconium alloys 
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eation in Al-Cu alloys [25l. |26|. 



29]. 



We should, though, note that here only the ordering (orientation) aspect of the issue 
is analyzed. This is characterized by a non-conserved order parameter field variable. The 
effect of the composition field governed by conserved kinetic equation is decoupled from 
the TDGL model and is not treated here. A more general set-up with coupled conserved 



and non-conserved field variables with defects in an elastic solid was presented in 



30|. The 



ground state exact solution with a mean field treatment of the static non-conserved field 



3l|. 



regarding second-phase nucleation on an edge dislocation was given in 

The organization of this paper is as follows. In Sec. [TT1 we lay down the formalism of 
the employed model, i.e. the governing TDGL equation and the parameters entering it for 
the problem under consideration. The numerical method utilized to solve the governing 
equation is outlined in Sec. IIHI The results of our numerical computations regarding the 
spatial-temporal evolution of the order parameter in the vicinity of an edge dislocation, for 
different values of the phenomenological coefficients in the TDLG equation, are presented 
in Sec. IIVI The steady-state solutions and the effect of the selected boundary conditions 
are also discussed in this section. The influence of the background thermal noise on the 
solutions is evaluated and discussed in Sec. [V] Section |VT] presents an analytic computation 
for late time evolution of the order parameter using Mazenko's approach. In Sec. IVIfl 
we give a summary of our main results, state our conclusions and also make some further 
remarks concerning generalizations of this work and applications to the real-life second-phase 
precipitation. 



II. MODEL DESCRIPTION 

We describe the kinetics of second-phase nucleation on an edge dislocation in a crystalline 
elastic solid by a time-dependent Ginzburg-Landau theory. A scalar order parameter field 
ri(p,t) depending on space p and time t defines the symmetry of the structure and dis- 
tinguishes the two prevailing phases, the solid solution and the precipitated second phase. 
Thus, t] = describes the high-temperature solid solution and t] ^ the second-phase 
nucleus. The non-conserved order parameter f](p,t) obeys the TDGL or Ginzburg-Landau- 



Khalatnikov (GLK) equation with a thermal noise i9 a 32|, |33(, which corresponds to model 



A in the classification scheme of Hohenberg and Halperin 34| . viz. 



dri _ 5 J 7 n . . 



Here, T a is a kinetic coefficient that characterizes the inter 



the total free energy of the system expressed by 
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35] 



(1) 

ace boundary mobility and T is 
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9a 

2 
K 



(Vr/) 2 + V{rf) + aV ■ urf 



+(t-t)( v ' u ' 



MVu:Vu 



dp. 



(2) 



Here g a (Vri) 2 arises from the inhomogeneity of rj and is related to short-range interactions, 
with g a taken as a positive constant. This term accounts for the existence of interfaces 



within an equilibrium inhomogeneous system 
Landau potential j^l in the form 



36|. The second term in the integrand is the 



V(?7) = ^r if + -u rj 4 + ^v Q r] e 



(3) 



where the coefficients r$ and uq are temperature-dependent variables while the term vor/ 6 
with v o > is needed for the stability when m < 0. The term arj 2 V ■ u in Eq. ([2]) describes 
the interaction between the order parameter and the displacement vector field u, where the 
strength of the interaction is designated by a and is assumed to be a positive constant. The 
second line in Eq. (EJ) is the elastic free energy, where K and M are the bulk and shear 
modulus, respectively; d is the space dimensionality; Vu = (duj/dxi + dui/dxj)/2 is the 
strain tensor, and the symbol : denotes the tensorial product contracted on two indices. The 
space integral in Eq. ([2]) is over the volume of the system. 

The function i9(p, t) in Eq. (pQ) stands for the presence of background random ther- 
mal motion, the so-called Langevin noise, prevailing at temperature T and satisfying the 
fluctuation-dissipation condition 

(0 a (p, t)0(p', 0) = 2k B TFJ(p - p')5(t - 0, (4) 
(0 a (p,*)) = O (5) 

where the averages (...) are over a Gaussian distribution function representing a Gaussian 
white noise, and ks is the Boltzmann constant. In this context, one quantity of interest in 
our study is the temporal evolution of the ensemble average of the square of the local order 
parameter, namely 

S(t) = ( V \p,t)), (6) 



where the average is taken over random initial conditions (noise) described by a Gaussian 
distribution function. 

Here, we assume that mechanical equilibrium for the displacement <5J r /5u = is satisfied 
at all times. A dislocation generates local strains that change the equilibrium condition in 
the solid. The mechanical equilibrium equation that includes the force field generated by an 
edge dislocation in the xy-p\ane with the Burgers vector in the x direction b = b e x is 

MV 2 u + (A - M) VV ■ u + a Vr/ 2 = -Mb e y 5(x)5(y) (7) 

where A = K + 2M(1 — 1/d), e y is the unit vector along the y-axis and £(•) is the Dirac 



delta 
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]. Equation (J7J) is then used to eliminate the elastic field from the expression 
for the total free energy (e.g. {mJ), which now can be expressed as 

f (Vr/) 2 + -r lV 2 + -u lV 4 + -v oV 6 ] dp. (8) 

The last three terms in the integrand correspond to the Landau potential given by Eq. (3) 
but with altered coefficients for the quadratic and quartic terms to account for the presence 
of an edge dislocation and an elastic body, respectively: 

r% = \r \ (sgn(r ) -p cos9/p), (9) 
U\ = uq — 2a 2 /A, (10) 

where po = 2abM/ (n \tq \ A) and b is the magnitude of the Burgers vector. Figure [1] shows 
the geometry of the dislocation. The parameter po can be considered as a local characteristic 
length related to the presence of the defect in elastic body. For a defect free crystal, r 1 = r , 
and for a rigid crystal, u\ = uq. In more detail, we may rewrite Eq. (Q in the form 

ri = a[T-T c (p,9)], (11) 
_ . „. _ 2abM cos 9 , . 

T ^^^K— ■ < 12 > 

Here a is taken to be a positive constant, T c q the phase transition temperature in a defect free 

crystal and T c (p, 9) the phase transition temperature for a crystal with an edge dislocation. 

The governing equation for the space-time variation of the order parameter is now obtained 

by inserting Eq. ([8]) into Eq. (JT|): 

L on 

fr~^ = 9aV 2 r] - (r lV + u^ 3 + v Q rf>) + o (p, t). (13) 

In this study, computations of phase transformation are performed with different combina- 
tions of r and U\ by solving Eq. ([TBI . 
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III. NUMERICAL METHOD 



We have used the open-source partial differential equation solver package FiPy [39] for 
our numerical computations. FiPy utilizes a standard finite volume approach, which is 
extensively used in computational fluid dynamics, in order to reduce the model equations to 
a form tractable to linear solvers . The spatio-temporal evolution of 77 in a two-dimensional 
space is computed, which is adequate, since a straight edge dislocation is considered here. 
The dislocation is placed at the origin, (x,y) = (0,0), with the slip direction along the line 
x = 0, see Fig. HJ We in general use a square mesh consisting of 200 x 200 equally-sized 
square elements, otherwise it is specified. Each element has a side length Al = po/20. 
The gradient of 77 perpendicular to the outer boundaries of the mesh is set equal to zero, 
i.e. n • V?7 = 0, where n is a unit vector perpendicular to a boundary, and thus periodic 
boundary conditions are achieved. In the cases where no thermal noise is included, the 
initial value of the order parameter, r]i nit , is taken to be a small positive random number 
of the order of 10~ 2 of the maximal value of 77 obtained in the computations. is a 

uniform rectangular distribution varying between 0.005 and 0.01. The system is considered 
to be large enough so that the periodic boundary condition would not affect the results. A 
reference time increment is defined as At re f = 0.9Al 2 /(2g a T a ). This time step is chosen to 
be sufficiently small to provide stable solutions for 77 for all the different combinations of tq 
and u\ studied here. 



IV. RESULTS 

The spatio-temporal evolution of 77 is presented for different sets of input parameters: 
{r , iti, v , g a }. First the situation when r > 0, i.e. T > T c0 , is investigated; for which, 
we set ro = 1. The elastic interaction energy, embedded in u\, is varied. Only results for 
cases with negative values of u\ are presented here; u\ = —1,-2,-3, respectively. This 
choice will be discussed in this section. The value of the coefficient of the Laplacian term 
in Eq. (fT3"l) is set as g a = 0.1, implying that the influence of the gradient of 77 is relatively 
large. The influence of g a will be also be discussed here. Thereafter, the behavior for r < 
is studied, namely, r = — 1 with u\ = {— 1, 1}. In all the calculations, v is put equal to 
unity, otherwise it is specified. It should be emphasized that the purpose of this study is to 
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explore the characteristics of the ordering evolution in the parameter space {ro,u±} G R. 



A. Spatio-temporal evolution 

Defects shift the nucleation temperature locally to a higher value and make it space- 
dependent; see Eq. (fl"2l) . Let us first consider the situation with T > T c0 , i.e. r > 0, 
and u\ < 0. Figure [2] shows the evolution of 77 when ro = 1 and u = —1. As can be seen, 
a peak emerges in the 77-surface, which is mainly situated on the compressive side of the 
dislocation (x > 0), and it evolves until it finds a stable shape, i.e. a steady-state solution 



is obtained. In Fig. 3(a), the temporal evolution of 77 at y — 0, i.e. for a cross-section along 
the x-axis, is shown at different times. A relatively sharp peak emerges first, and afterward 
some broadening occurs. The maximum value of this peak, r] max = max [77 pea A;], is reached at 
t pd 1200At re /, and the steady-state shape comes after at t ~ 1800Ai re j. It is observed that 
after a relatively slow start, the evolution rate increases until the steady-state value, r] max , is 
nearly achieved. Afterward, it decreases just before the system reaches steady state. Next, 
the case with u± = —2, tq = 1 is investigated, and the evolutions of the order parameter are 



presented in Figs. H] and 3(b) It is seen that a peak first develops reaching its maximum 
value, rj max , at t ps 1200At re j, as in the case with ui = —1. Then a considerable broadening 
takes place until a steady-state shape is reached at t m 7000 At re f. We have also studied 
the situation where u\ = —3 and tq = 1. The 77-surface at three different times displayed 



in Fig. and Fig. 3(c) shows r](x,y = 0) for every 200 time-step. The order parameter 
evolution here differs from the cases with u\ = —1 and u\ = —2. Again a peak develops that 
afterward broadens, but now the second phase grows until the whole material is transformed 
to this phase except near the dislocation on its tensile side (x < 0) where the transformation 
is held back. The peak reaches its maximum value, r] max , at t ~ 1200At re j, and steady state 
is obtained at t ~ 5000 At re f. A close-up of the 77 field, in the form of a contour plot, is 



shown in Fig. 5(d) 



The case r = — 1 (T < T c0 ) corresponds to a quick decrease of temperature, i.e. a 
quench below the transition temperature of the undisturbed system. At the same time as the 
temperature changes, an edge dislocation is introduced. First, the results from computations 
with Ui = 1 are presented. As can be seen clearly from Fig. [61 a peak is growing where 
compressive stresses are induced by the dislocation (i.e. x > 0) with its maximum close to 
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the dislocation. The evolution of rj is also shown in Fig. 3(d) as profiles of f](x,y = 0). 
Here, the peak evolves until t « 300At re j. Subsequently, t] increases in the entire material, 
while leaving a valley on the tensile side of the dislocation. The time to reach steady state 
is found to be t ~ 900At re j. 



For ro = — 1 and u\ = —1, the development of rj has the same character as with u\ — 1, cf. 



Fig. 3(e) The steady-state value of rj is, however, larger in the whole plane, while maximum 
values are reached at approximately the same time as with u\ — 1. Steady state seems to 
be achieved at t m 900At re j. The steady-state solutions for ro = — 1 is similar to that of 
ro = 1 and u\ = —3. However, nucleation of second phase would not initiate in the tensile 
region of the solid. Instead, for r < 0, after the local phase transition near the dislocation, 
the order parameter rj increases in the whole material at the same time. 



To recap, the values of some characteristic measures for the studied cases are given in 
Table [H It is seen that r] max increases with decreasing ui for both positive and negative values 
of ro- Equation UM tells us that u\ decreases by increasing the interaction strength a, i.e. the 
interaction between 77 and u. Thus, with a larger value of a, the singularity of the stress field 
induced by the dislocation, contributes to a larger value of rj. Also, decreasing the elastic 
moduli, i.e. decreasing A, reduces u\. The time needed for r] pea k to reach its maximum rj max 
is denoted by t mp . In the case of r = 1, a top that emerges near the dislocation is found 
to grow four times slower than in the case of quenching, i.e. with r = —1. This could be 
explained by indicating that in the quenching situation a defect-free material is expected 
to fully transform into the second phase, and thus the driving force for phase transition is 
greater for a system with r < 0. Correspondingly, the time to reach steady state, t ss , is 
shorter with ro = — 1 than with ro = 1 for u\ = —1. For the latter case, the evolution of 
a relatively thin top is faster than for lower values of U\. With u\ = —2, a larger area at 
the dislocation transforms and it takes about four times longer. When practically the whole 
material is transformed into the second phase, u\ = —3, the growth rate increases, which 
may be expected. To estimate 77 as p — > 00, the value at (x,y) = (5poj5po) is used as a 
measure and is denoted by r/^. The results are discussed further in subsection IIV CI 
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B. Nucleus shape evolution 



To illustrate the evolution of a second phase in the (x, y)-plane, black and white patterns 
representing the two phases are displayed in Fig. [7J Black indicates areas where r\ > r] max /2 
and represents the second phase and vice verse. This choice is somewhat arbitrary, but 
renders patterns that illustrate the evolution of the shape of a nucleus. Patterns for four 
different combinations of r and U\ at three different times are given in the figure. In the 
upper row, Figs. 7(a) to 7(c) , patterns for the case r = 1 and U\ = —2 are shown. It 



is clearly seen that the second phase evolves in the shape of a circle and is located on the 
compressive side of the dislocation, i.e. where x > 0. For U\ = —3, the growing second phase 

in 



is not circular as in the other cases. Instead, it is kidney shaped, see Fig. 7(d) and 7(e^ 



the second row. As steady state is reached, see Fig. 7(e) , the area where the transformation 
is held back by the dislocation is circular. 

When r = —1, the patterns for U\ = 1 and U\ = —1 are shown in the two lower rows of 
Fig. [7J respectively. It is seen that the second-phase evolution is somewhat slower for x < 0. 
The only significant difference in the behavior is the way that the matrix phase shrinks; cf. 



Figs. 7(g) to 7(i) with 7(j) to 7(1) In steady state, the area where the transformation is held 
back is circular in both cases, although it is smaller with the negative u%. 



C. Steady-state solutions 

Let us consider the long-time limit when steady state is reached, i.e. drj/dt = 0. Fur- 
thermore, we assume the system is away from the critical region, where the long wavelength 



TABLE I. Comparison between results {rj max , t mp , t ss , r/oo} for different input values {ro, u{\. 



r => 


1 


1 


1 


-1 


-1 


U\ =>• 


-1 


-2 


-3 


1 


-1 


1]max 


1.46 


1.78 


2.05 


1.37 


1.77 


tmp 1 ' ^t re j 


1200 


1200 


1200 


300 


300 


tss 1 ' At re f 


1800 


7000 


5000 


900 


900 


r]oo 








1.63 


0.81 


1.29 
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fluctuation in the order parameter can be neglected, i.e., g a V 2 f] ~ 0. Then ignoring the 
noise term, Eq. f fT3|) gives: 

f\ (r x + uiff + v fj 4 ) = 0, (14) 

where f] is the elastic medium (mean-field) order parameter containing the edge dislocation. 
This reduced equation corresponds to dV{ff)/dfj = with modified coefficients r 1 and Ui 
instead of r and u , cf. Eqs. ([8]) to ([TO]) . The non-zero solutions of Eq. (fl~4l) are 

fjl = - Ml± ^- 4 "° ri . (15) 

The solutions f] + exist for r\ < uf/4vo, whereas those of f)_ exist for < r\ < u\/Avq. We 
also note that V{ff) = together with dV(fj)/df] = yield r\ = 3uf/lQvo. 

Figure M illustrates two different settings of the the coefficients of V{rf): (a) U\ > 0, and 
(b) U\ < 0, both with v > 0. If u\ > and r\ > 0, no nucleation of second phase would 
occur, whereas ri = sets off the nucleation and r\ < corresponds to formation and an 
unstable growth state of the second phase. In the case of u\ < 0, at high temperatures, 
ri is large and positive and V(ff) has a simple structure with a minimum at 77 = 0. As 
the temperature is decreased r\ gets smaller and in the range r\ < u1/4vq, u± < spawns 
two local minima at f] 7^ 0, i.e. emergence of metastable second-phase states, while the 
solid solution is stable (the global minimum). A further decrease of temperature makes 
ri < 3M^/16t>o, at which the global minimum shifts from 77 = to the two symmetrically 
located states with f] ^ 0. This marks a first-order phase transition. For r± < the local 
stability of the solution fj = disappears. 

The parameter r\ is a function of both temperature and spatial coordinates. Accordingly, 
so are f] + and the location where nucleation and further evolution of the second phase would 
occur. Depending on the sign of u\, different situations would arise. Figure H] shows 77 along 
the x-axis obtained by the numerical computations, r] num , which includes the term g a V 2 ?7 
with g a = 0.1, and the aforementioned analytical solution, r] ana = fj, for different parameter 
settings. We see that the main features of t] ana can also be found from the numerical 
solutions. The deviations are largest near the dislocation line, where r] ana goes to infinity, 
and at the locations where there are sharp boundaries between phases. This is expected 
since the phase-field model (numerical solution) produces smooth interfaces. For situations 
with u\ > 0, regardless of magnitude, Eq. ()9]) tells us that the ordered phase would evolve 
within a circle confined by p = p cos9, if r > 0. In the quenched case (r < 0), the whole 
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material would transform except where p < — po cos 9, i.e. inside a circle of a radius po- 



The values of u± would only influence the magnitude of r] + . The profiles r\ ana in Figs. 9(a) 



and 9(e) show that boundaries between phases in both cases are continuous (second-order 
transition). For negative ux, the radii of the limiting circles vary with ui as Po/(l — 3n^/16), 
since a first-order transition occurs if r\ < "5u\/\Qvq. In Figs. 9(b) to 9(d) and 9(f), it is 



seen that the order parameter is discontinuous at the corresponding locations. 

We should note that the smoothness of the interfaces in the numerical computations 
is governed by the term g a V 2 ?7 in Eq. ( ITHj) . For instance, if the value of the coefficient 
g a is reduced by two orders of magnitude, the softness of the phase interface is reduced 
substantially, giving rise to a sharper shape that is closer to rj ana . 



D. Influence of periodic boundary conditions 



In the case of r = 1 and U\ = —3, cf. Fig. 5(c) and 9(d), the periodic boundary 
conditions inevitably have an impact on the second-phase evolution; the broadening of the 
top continues until the boundary is reached at t ~ 1800At re f. This means that there is 
an interaction with the corresponding dislocations that are mirrored at the boundaries, i.e. 
those located at (x, y) = ±10po However, if the mesh size is increased, thus increasing the 
distance to the image dislocations, the same behavior is obtained. That is, no significant 
difference in growth rates is found, except that the time that would take for the second 
phase to come into contact with the edge of the grid would be twice as long. The relatively 
large difference between r] ana and r) num for the case tq = 1 and u\ = —2 is probably due to 
the imposed periodic boundary conditions, but that impact was not investigated further. 



V. INFLUENCE OF THERMAL NOISE 

In the preceding computations, the Langevin thermal noise, i.e. the d a term on the right 
hand side in Eq. (1131) . was disregarded. In order to investigate the influence of this term 
on the temporal evolution of the system, numerical computations have been performed for 
situations with and without the presence of a dislocation. The thermal noise is represented 
by a Gaussian distribution with the mean value equal to zero and the variance of 10~ 4 . 
Instead of using a uniform random distribution as the initial value for rj, cf. in section IHfl 
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a constant initial value equal to 7.5 x 10~ 3 is used. Then for each time increment a new 
noise distribution is added. Here, we have used a finer mesh, namely a mesh with 400 x 400 
elements, compared to our foregoing computations, and also the time increment is taken to 
be smaller than before, i.e., only a fourth of the previous one: At = At re f/A. 



A. System without a dislocation 

For reference, a system without a dislocation and U\ = — 1 is studied. The system is 
quenched from above to below the critical temperature T c instantaneously, i.e. from tq > 
to ro < 0. In such a case the whole system is expected to transform into the second phase, 



see e.g. Fig. 8(b) First a computation without the thermal noise term is performed. The 
top row of Fig. [10] shows the evolving domain at three different periods. It is clearly seen 
that the domain is coarsening with time. We have not yet quantified the level of coarseness 
as time advances. It should be mentioned that the white and black areas do not represent 
the matrix phase and second phase, respectively, instead a mean value of r\ over the whole 
area, r] mean , at each time step is determined, and locations with rj greater than this mean 
value are indicated in black, and viceversa. Since the GLK equation is non- conservative, 
the entire system is eventually engulfed by the second phase. From the plot of r] mean versus 
t the variation in the rate of transformation can be inferred, Fig. [TTJ It is seen that after 
an initial period with quite a slow increase in t], the transformation accelerates until the 
equilibrium value is reached. The curves in Fig. [TTJ fit a logistic curve of the form 

V(t) ~ - t^—, n, (16) 

/w l + 6exp(-ct)' v ; 

with a = 1.289, b = 5486.75, and c = 0.02242. The fit is quite good for t/At ref > 400. 

The corresponding calculation is done with the presence of thermal noise as described 
above. In Fig. [10] the bottom row shows the obtained patterns for the same times as for 
the case without noise. A somewhat coarser domain is obtained for the situation with 
thermal noise. No significant difference in the rate of phase transformation is, however, 
found, see Fig. [TTJ Based on these observations, we conclude that the influence of thermal 
noise is insignificant for the domain growth on quenching in the absence of dislocation, at 
least for the parameter set used in our comparative study. We shall recall, though, that 
the absence of the thermal noise term in Eq. f[T3l assumes that the system is effectively 
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evolving at T = K. In this context, it has been pointed out that the renormalization 
group (RG) theory shows three fixed points for temperature, namely, T = 0, T = T c and 
T = oo 40j. This implies that during ordering, temperatures below T c move toward zero, 
while those under disordering move toward infinity in RG flow maps. Therefore, the final 
quenching temperature Tp is an impertinent variable for quenches into the ordered phase, 
which justifies the neglect of thermal noise in the equations of motion j4o| . 

The domain shapes observed in Fig. are equivalent to the snapshots of coarsening 



produced by mode 
KIM model 
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s with non-conserved order parameter dynamics, e.g. the kinetic Ising or 
4l| , which is in the same university class as the GLK model employed here. 
In more detail, the recent publication by Olejarz and colleagues ^1], who used a 2D KIM 
model on a 1024x1024 square lattice with periodic boundary conditions, show very similar 
snapshots of coarsening as our Fig. [TT], after a quench from T = oo to T = 0. This again 
vindicates our conclusion regarding the effect of thermal noise on the kinetics of domain 
growth. The time evolution patterns in Fig. [TUlmay also be compared with those produced 



in Onuki's book 



42| , where the results of an exact numerical solution of model A (with cubic 



order parameter) are compared with other known approximate solutions in the literature. 



B. System with a dislocation 

We now study the influence of the thermal noise term in Eq. (fT3]) on the phase transfor- 
mation kinetics in the presence of a dislocation for one case; u\ = —1 and tq = 1, for which 
results of computations without the noise term were presented earlier, see subsections IIV Al 
Here, the evolution of the order parameter is illustrated in the form of contour plots of rj 
at different times, see Fig. [121 It is seen that at a very early stage of phase transition a 
top emerges near the origin in the half-plane where x > 0, although the variance of the 
background noise is of the same order of magnitude. With increasing time the thermal noise 
loses its impact. A single top is expected to evolve, which is confined within a circular area 
on the compressive side of the dislocation, as in the case with no noise added, cf. Fig. [21 
However, as can be seen from Fig. [THJ, there is a difference in the evolution rate if noise 
is included in the calculations. For the disturbed system, f] pe ak evolves faster, in the sense 
that the growth of the top is larger in the very beginning of the evolution. At later times, 
the growth rate is very similar between the two cases as can be seen from Fig. [131 The 
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two curves will coincide if the dashed line is shifted horizontally by about 300At re j. But if 
the variance of the Gausssian distribution is reduced, this shift is also reduced. We should 
also note that the nucleation rate of the second phase on dislocation is much higher than in 
homogeneous (defect free solid) nucleation. This is because the nucleation barrier energy at 
an edge dislocation is lower than that for homogeneous nucleation; see e.g. jsj]. A similar 
remark also applies to the growth of second phase 43|. 

It can be concluded that the results in this study will not differ much if Langevin thermal 
noise is included in the governing equation or not. The only significant difference that 
has been noted is that the nucleation of the ordered phase at a dislocation starts earlier if 
thermal noise is taken into account at T > T c q, i.e. an edge dislocation is hosted in the solid 
solution. This can be interpreted as that the noise supports the nucleation at a singularity, 
and may not suppress it. The main features of the evolution of all the studied cases in this 
paper will still be effective in the presence of thermal noise, at least with the parameter set 
and noise distribution used here. 



VI. ANALYTICAL SOLUTIONS: LATE TIME GROWTH 

It is worthwhile to discuss now the late time growth behavior of the second-phase/matrix 
interface within the GLK theory. To this end, we consider the equal-time order parameter 
correlation function 

C(p,t) = ( V (p,t) V (0,t)). (17) 
' n late times after a deep temperature quench into an unstable state, C(p,t) takes a form 

C(p,t)= V 2 F(p/L(t)), (18) 

where rjo = rj + is the magnitude of long-time equilibrium value of rj(p,t) and F(») is a 
universal scaling function, which depends on the symmetry of the system (here a scalar 
order parameter) and the spatial dimension. Here, L(t) is a characteristic length scale 
featuring the growth law for the second phase. 



Mazenko's functional-integral method for phase ordering kinetics 45|, |46] shows that the 
scaling function F(») obeys the eigenvalue problem 

- /ix ■ V X F = tan(7rF/2) + V 2 X F, (19) 
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TABLE II. Computational results of late time temporal evolution of the mean order parameter 
using Eqs. (J2T]) with up to n = 3, and ([H with d = 2, cf. Fig. [HJ 



T a t VS/ m (A = 0.1) VS/ Vo (A = 1) 

10 0.983266 0.912301 

100 0.994734 0.972326 

1000 0.998337 0.991307 

10000 0.999475 0.997258 



where x = p/L(t) and [i is an eigenvalue to be determined. An exact analytical expression 
for L(t) valid for all times is not known, however, at late times 

L(t) = (AT af it) 1/2 . (20) 



en 48] curvature-driven growth 



46j has determined numerically 



Equation (120]) is the manifest Lifshitz 47|, Cahn and Al 
law for non-conserved systems at late times. Mazenko 
the eigenvalue /i = fi*(d), which is space dimension dependent, namely, /i = 1.104 and 
fi = 0.5917 for d = 2 and d = 3, respectively. 

Let us now using Mazenko's method calculate the temporal evolution of the ensemble 
average of the order parameter (r] 2 (p, t)) 1 ^ 2 , cf. Eq. ([6]), at the interface for r\ < and 



Mi < 0. The computations outlined in Appendix lAl give 



n=0 

where a^'s are constants for a given state of the system. This series is asymptotically 
convergent. In Appendix [Al we have given the first three coefficients of the series fl2"T|) . which 
are adequate for convergence. Figure [TH depicts \fS /r] as a function of the dimensionless 
time for two values of A = \ri\v /u\ using Eqs. (121]) with up to n = 3, and (|20!) with d = 2. 
It is seen that an increase in A slows down the completion of equilibrium. Table [TT] shows 
this quantitatively. 

It may be of interest to compare the results of the late time evolution of (r] 2 (p, t)) 1 ^ 2 
with the numerical computations of r] mean made on the entire range of t depicted in Fig. 
[TT1 and discussed in Sec. IV Al for dislocation-free solid. As mentioned there, the curves in 
Fig. [TT] fit the logistic curve, given by Eq. f lT6|) . fairly well. Although, we do not know 

1(3 



where exactly to draw the line for late time growth and even the two quantities may not 
directly be comparable, we have attempted to fit the calculation output for A = 1, Fig. [TH 
to a logistic curve, Eq. ( I16p . in the time interval T a t G [0.5, 10]. The fit was adequate with 
a = 0.902,6 = 0.386, c = 0.386, and t -»■ T a t for S 1 ' 2 /^ = (r] 2 {p, t)) 1/2 /r/ - This shows that 
our numerical solutions for r] mean follow the same evolution trend as the analytical solution 
using Mazenko's theory. 

VII. SUMMARY AND CONCLUSIONS 

In this work, we have used the time- dependent Ginzburg-Landau (TDGL) equation for 
a single component non- conservative structural order parameter rj to model the spatio- 
temporal evolution of a second phase in the vicinity of an edge dislocation in an elastic 
crystalline solid. A symmetric Landau potential of a sixth-order, r/ 6 , was employed. The 
phase-field equation was solved numerically using a finite volume method, where a wide 
range of parameter sets is explored. Computations were performed for the situations where 
the temperature was held above the transition temperature of a defect-free crystal, T c q, as- 
well- as below it. In both cases, the influence of the elastic properties of the material and 
the strength of interaction between the order parameter and the elastic displacement field 
were examined by varying a model parameter that comprises these properties. 

We found that the introduction of a dislocation always triggers nucleation of a second 
phase. The phase transition initiates in the vicinity of the dislocation line in the region 
where stresses are compressive, regardless of the parameter setting. If the temperature is 
above T c0 and the elastic interaction is moderate, we found that the second phase grows 
locally within a confined space. However, if the elasticity and/or elastic interaction is large, 
eventually the phase transition will spread throughout the whole material even though the 
temperature exceeds T c0 . In the regions where the stresses induced by the dislocation are 
tensile in character, the phase transition is suppressed, though not fully. If the system 
is quenched below T c o, the entire material transforms and steady state is first reached on 
the compressive side of the dislocation. On the side with tensile stresses, the evolution of 
ordering is held back. The steady-state distribution of phases is estimated by considering a 
modified Landau type potential. 

The influence of the Langevin thermal noise term in the TDGL equation was also exam- 
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ined. We found that if the dislocation is introduced in the crystal above T c o, the thermal 
noise will support the nucleation of the ordered second phase, and steady state will be 
reached earlier than if the thermal noise were absent. But if a dislocation is introduced in a 
solid solution whilst the system is being quenched below T c0 , the evolution of second phase 
ordering will not be affected significantly by the noise term. 

For a dislocation-free solid, we compared our numerical computations for a mean-field 
(spatially-averaged) order parameter as a function of time with the late time growth of the 
ensemble- averaged order parameter calculated from Mazenko's theory of domain growth, 
and found that both results follow late time logistic curves. 

The present work is part of a larger study on how singular stress fields affect phase 
transformations in crystalline materials, where both dislocations and cracks are considered. 
The subsequent steps of this study aim to extend the present calculations to include a two- 
component field structural order parameter, and to consider the coupling of the composition, 
obeying a conserved kinetic equation, and the two-component order parameter. This delib- 
eration not only will describe the formation and growth of second phase in the vicinity of 
defects but also will tell its orientation in the presence of applies stress. 
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Appendix A: Evolution of the interface in quenched phase 



Here we outline the method introduced by Mazenko and coworkers |45l. l46l. |49| to study the 
growth kinetics of quenched systems from early through late times which comprised sharp 
interfaces in the late stage development. We tailor the approach to our application, i.e. 
the growth of a second phase within the framework of non-conserved rf Ginzburg-Landau 
theory. 

We consider the time-dependent Ginzburg-Landau equation for an elastic body containing 
a dislocation. As in Sec. |V] the considered equation is supplied by an initial probability 
distribution governing 77 at time t = in the solid solution T > T c , which is supposed to 
be Gaussian with the initial correlation (rj(p, 0)r/(p', 0)) = 2/csT5(/o — p') and (t](p, 0)) = 0. 
The system is then quenched to a state below T c upon which a second phase is formed 
around the dislocation. One can assume that the system is quenched to zero temperature 
so the thermal noise can be set equal to zero for T <T C . 

A key notion in Mazenko's theory is the separation of the order parameter field 77 into a 
peak contribution a and a fluctuating term drj, viz. 

r)(p,t) = a[s(p,t)} + 6r](p,t), (Al) 

where a[s(p,t)} is a functional of the auxiliary field s(p,t) and s is given a physical in- 
terpretation, near interfaces (walls), as a coordinate normal to the wall. So a[s] may be 
considered as the equilibrium interfacial profile depending on s. Furthermore, it is assumed 
that s is proportional to the characteristic length L that scales the system. More precisely, 
one defines 

5 2 = (s 2 (p,t)) = L 2 /n, (A2) 

where the angular brackets denote the ensemble average over random initial conditions, and 
also one posits that s and L are increasing function of time. 

Ignoring now the fluctuating term in Eq. flAlj) , then placing it into the expression for the 
free energy functional, Eq. (j5J, and using the Euler-Lagrange equation, we obtain 

|^-(-hk + «/+V 5 )=0, (A3) 

which is precisely the time-independent Ginzburg-Landau equation near an isolated defect, 
cf. Eq. ffTB"]) . We can rewrite this last equation in terms of scaled (dimensionless) variables 
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in the form 



ld 2 a 



2 ds 2 



a - sgn(ni)o- 3 - Act 5 = 0, 



(A4) 



where a = a/a*, a* = (Iril/liiil) 1 ' 2 , s = (ri\/ g a ) l l 2 s, and A = \ri\v /u\. We now, for 
convenience, drop the bar and the hat symbols from the variables and rewrite Eq. ( 1A4I) in 
the form 

ldV 



V'[a] = -a + sgn(wi)cr 3 + Act 5 . 



Moreover, we impose the boundary condition 



hm — = 0. 

s->oo ds 



Integration of ( 1A5|) gives 



dx 



2 io y/V[x]-V[ Vo ] 



where the critical points ir/ 2 are found from V'[?7o] = 0, 



V 2 o 



^-(-sgn(u 1 ) + v / l + 4A). 



Furthermore, 



We evaluate the integral in (|A8[) to obtain 

s = ±£o ^ arctan 6 + i— j , 
with Co = (2 - sgn(n 1 )r/o)" 1/2 > 

2 = ^o 2 (% 2 + l + ^(% 2 + 4)). 

Inverting Eq. flAllj) to find cr as function of s, after some manipulation, we obtain 



a 



77o\/l — etanhs 
\/l — e tanh 2 s 



2 /^ 2 + 1\ 

with e =3(^J for u 1 <0, 

S = s/6)- 



(A5) 
(A6) 



(A7) 
(A8) 

(A9) 
(A10) 

(All) 
(A12) 
(A13) 

(A14) 

(A15) 
(A16) 



Equation (1A14I) is identical to Mazenko's result 46[ for uq > except for some modifications 
in the definitions of e and £ - 



22 



It can be shown that as A is increased a tends toward the origin. Moreover, for large s, 
a ~ tjq. This trend can seen by expanding V[x] near rj and evaluating the integral in Eq. 
flA"8l) to obtain 

a[s] = ±77o (l - exp[-V2V"(77o) «])■ (A17) 

We should also recall that for «i > the case A = gives rjo = ±1 with a = ±tanh[s]. 
One quantity of interest for our analysis is the interfacial width £ (46| defined by 



i r-t-oo 

t=S K-^NJds. (A18) 

Substituting for <j(x) in the integrand from Eq. (IA14I) and evaluating the integral gives 

£ = 2^=arctanh[v^]. (A19) 
v e 

Using a standard identity in hyperbolic functions, we write this formula as 

£ = log 1 + ^1 , for < e < 1. (A20) 
V e 1 - V e 

Let us next calculate the evolution of the local order parameter, cf. Eq. (jSJ), in the s- 
coordinate 

S(t) = (a 2 [s}), (A21) 

with 

/+oo 
(j 2 [x}V(x)dx, (A22) 
-00 

where V(x) is taken to be a Gaussian distribution of a real random variable x with a variance 
s, defined in Eq. (1A2I) . given by 

P(x) = -=L=exp[-x 2 /2s 2 ], (A23) 

V27TSf 



Again following the procedure delineated in 



46], we split 



a 2 [s] = a 2 [oo] + Aa 2 [s], (A24) 
where A<r 2 [s] = & 2 [s} — rj 2 .. Hence Eq. (1A21|) is expressed in the form 



i r+oo 

S(t) = a 2 [oo] + -== \ e- x2/2s2 Aa 2 [x}dx, (A25) 

V27TS 2 J-oo 

Expanding the exponential term in Eq. (IA25j) in infinite series, 



1 "1 i /*H~oo 
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where we put <7 2 [oo] = rfo by (1A17I) . Considering the first three terms in the expansion, then 
performing the integrations and replacing s by L(t) via Eq. (IA2[) . we obtain 



s W = „2(i-^ + ||-§ + ... 



where a± — £/ a/2, and 



71-6 Zfrt , 71 & 



2c2 




€ 



a 3 = ^|(r + -^), (A28) 



«5 



2 5 / 2 80 



(e 4 + ^ 2 + ^). (A29) 



The series (IA270 is asymptotically convergent; the expansion up to 0(L~ 5 ) provides sufficient 
accuracy in computations. Hence, if knowing the time- dependence of the characteristic 
length L, then the time evolution of the averaged order parameter can be calcuated from 
Eq. (|A27]). 

50| | has provided an alternative approach to compute the temporal evolution of the 



Cardy 



non-conserved order parameter in the case of a quench to a temperature at, or just above, the 
critical temperature, when the slowness in the dynamics is a result of the critical slowing 
down of local fluctuations. He found, using the renormalization group and e-expansion 
technique, to all orders in e, that the local fluctuations in the order parameter scale like 
t" 1 / 2 , and have a universal distribution. 
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Appendix B: Figures 




FIG. 1. Geometry of the edge dislocation. 
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yi p - 5 -5 xi p o yi p -5 -5 x , P() 

(c) t = 1000At re/ (d) t = 2000At re/ 

FIG. 2. Spatial distribution of r/ at various times t = 600,800,1000 and 2000 At re f, for ro = 1, 
u± = — 1 and vq = 1. At re f is the reference time-step used in computations and is denned in Sec. 
IIIII of the main text. Only every fifth node in the mesh is selected for the illustration. 
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(a) ro = 1 and u± = — 1 



(b) ro = 1 and iti 
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(c) ro = 1 and iti = 
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(d) ro = — 1 and u\ = 1 



-5 -4 -3 -2 -1 1 2 3 4 5 



(e) ro = — 1 and u\ = — 1 



FIG. 3. Evolution of the order parameter rj(x/po,y = 0) for different combinations of ro and 

with vq = 1. The time advances vertically in the plots; in a) to c) rj is given from t = 200At re f 

10000 At ref , and in d) and e) from t = 50At re/ ->■ 2000At re/ . 
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y/ P -5-5 x l p Q yl p Q -5-5 x / p 

(a) t = 1000At re/ (b) t = 2000At re/ 








(c) t = 10000Ai re/ 

FIG. 4. Spatial distribution of rj at various times t = 1000, 2000 and 10000 At re f, for ro = 1, 
ui = —2, and Uq- Only every fifth node in the mesh is selected for the illustration. 
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yi P " 5 "5 xl p Q 

(c) t = 6000At re/ (d) Contours at t = 6000At re/ 

FIG. 5. Spatial distribution of r/ at various times t = 1000, 2000 and 6000 At re f, for tq = 1, 
u\ = — 3 and vq. Only every fifth node in the mesh is seleted for the illustration. 
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(c) t = 2000At re/ (d) Contours at t = 2000Ai re/ 

FIG. 6. Spatial distribution of rj at various times t =200, 400 and 2000 At re f, for ro = —1, u\ = 1 
and vq. Only every fifth node in the mesh is selected for the illustration. 
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FIG. 7. Snapshots of evolution of a second phase (black) from the matrix material (white) for the 
cases (a)-(c) r = 1, u\ = -2, (d)-(f) r = 1, u± = -3, (g)-(i) r = -1, u\ = 1, and (j)-(l) r = -1, 
iti = —1. 
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(a) ui > (b) mi < 

FIG. 8. The Landau potential V{rf) = \r\if + \u\if + \v^rf with vq > versus 77 for various 
combinations of the coefficients. 
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FIG. 9. Comparison between analytical and numerical computations for spatial variation (x, y = 0) 
of order parameter in steady state for different parameter sets (ro, u±) with vq = 1. 
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FIG. 10. Snapshots of domain coarsening on quenching for the case of no dislocation. Top panel: 
without thermal noise term; bottom: with the thermal noise term. Here, tq = — 1, u\ = —1,Vq = 1. 
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FIG. 11. Temporal evolution of the system without a dislocation: mean field order parameter 
Vmean versus scaled time t/At re f with and without the presence of thermal noise (ro = —l,u± = 

-Mo = i)- 
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FIG. 12. Contours of rj when thermal noise is present for a system with a dislocation for the case 
ro = 1, U\ = — 1 and Vq = 1 at different times. 
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FIG. 13. Comparison of the temporal evolution of f] pea k with and without the presence of thermal 
noise for rg = 1, U\ = — 1 and vq = 1. 
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FIG. 14. Late temporal evolution of the root mean square of the order parameter scaled with the 
equilibrium value rjo at two values of the coefficient A. 
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